Shock Formation in a Multidimensional Viscoelastic Diffusive System

We examine a model for non-Fickian "sorption overshoot" behavior in diffusive polymer-penetrant systems. The equations of motion proposed by Cohen and White [SIAM J. Appl. Math., 51 (1991), pp. 472–483] are solved for two-dimensional problems using matched asymptotic expansions. The phenomenon of shock formation predicted by the model is examined and contrasted with similar behavior in classical reaction-diffusion systems. Mass uptake curves produced by the model are examined and shown to compare favorably with experimental observations

A vicinal surface model for epitaxial growth with logarithmic free energy

We study a continuum model for solid films that arises from the modeling of one-dimensional step flows on a vicinal surface in the attachment-detachment-limited regime. The resulting nonlinear partial differential equation, $u_t = -u^2(u^3+\alpha u)_{hhhh}$, gives the evolution for the surface slope $u$ as a function of the local height $h$ in a monotone step train. Subject to periodic boundary conditions and positive initial conditions, we prove the existence, uniqueness and positivity of global strong solutions to this PDE using two Lyapunov energy functions. The long time behavior of $u$ converging to a constant that only depends on the initial data is also investigated both analytically and numerically.Comment: 18 pages, 7 figure

Stability of Self-similar Solutions for Van der Waals Driven Thin Film Rupture

Recent studies of pinch-off of filaments and rupture in thin films have found infinite sets of first-type similarity solutions. Of these, the dynamically stable similarity solutions produce observable rupture behavior as localized, finite-time singularities in the models of the flow. In this letter we describe a systematic technique for calculating such solutions and determining their linear stability. For the problem of axisymmetric van der Waals driven rupture (recently studied by Zhang and Lister), we identify the unique stable similarity solution for point rupture of a thin film and an alternative mode of singularity formation corresponding to annular “ring rupture.

Stability and Dynamics of Self-Similarity in Evolution Equations

A methodology for studying the linear stability of self-similar solutions is discussed. These fundamental ideas are illustrated on three prototype problems: a simple ODE with finite-time blow-up, a second-order semi-linear heat equation with infinite-time spreading solutions, and the fourth-order Sivashinsky equation with finite-time self-similar blow-up. These examples are used to show that self-similar dynamics can be studied using many of the ideas arising in the study of dynamical systems. In particular, the use of dimensional analysis to derive scaling invariant similarity variables is discussed, as well as the role of symmetries in the context of stability of self-similar dynamics. The spectrum of the linear stability problem determines the rate at which the solution will approach a self-similar profile. For blow-up solutions it is demonstrated that the symmetries give rise to positive eigenvalues associated with the symmetries, and it is shown how this stability analysis can identify a unique stable (and observable) attracting solution from a countable infinity of similarity solutions

Gravity-driven Thin Liquid Films with Insoluble Surfactant: Smooth Traveling Waves

The flow of a thin layer of fluid down an inclined plane is modified by the presence of insoluble surfactant. For any finite surfactant mass, traveling waves are constructed for a system of lubrication equations describing the evolution of the free-surface fluid height and the surfactant concentration. The one-parameter family of solutions is investigated using perturbation theory with three small parameters: the coefficient of surface tension, the surfactant diffusivity, and the coefficient of the gravity-driven diffusive spreading of the fluid. When all three parameters are zero, the nonlinear PDE system is hyperbolic/degenerateparabolic, and admits traveling wave solutions in which the free-surface height is piecewise constant, and the surfactant concentration is piecewise linear and continuous. The jumps and corners in the traveling waves are regularized when the small parameters are nonzero; their structure is revealed through a combination of analysis and numerical simulation

Blowup and Dissipation in a Critical-Case Unstable Thin Film Equation

We study the dynamics of dissipation and blow-up in a critical-case unstable thin film equation. The governing equation is a nonlinear fourth-order degenerate parabolic PDE derived from a generalized model for lubrication flows of thin viscous fluid layers on solid surfaces. There is a critical mass for blow-up and a rich set of dynamics including families of similarity solutions for finite-time blow-up and infinite-time spreading. The structure and stability of the steady-states and the compactly-supported similarity solutions is studied

Stability of Traveling Waves in Thin Liquid Films Driven by Gravity and Surfactant

A thin layer of fluid flowing down a solid planar surface has a free surface height described by a nonlinear PDE derived via the lubrication approximation from the Navier Stokes equations. For thin films, surface tension plays an important role both in providing a significant driving force and in smoothing the free surface. Surfactant molecules on the free surface tend to reduce surface tension, setting up gradients that modify the shape of the free surface. In earlier work [12, 13J a traveling wave was found in which the free surface undergoes three sharp transitions, or internal layers, and the surfactant is distributed over a bounded region. This triple-step traveling wave satisfies a system of PDE, a hyperbolic conservation law for the free surface height, and a degenerate parabolic equation describing the surfactant distribution. As such, the traveling wave is overcornpressive. An examination of the linearized equations indicates the direction and growth rates of one-dimensional waves generated by small perturbations in various parts of the wave. Numerical simulations of the nonlinear equations offer further evidence of stability to one-dimensional perturbations